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Abstract 

Q H | We define a model of quantum computation with local fermionic modes (LFMs) — sites 

which can be either empty or occupied by a fermion. With the standard correspondence 
between the Foch space of m LFMs and the Hilbert space of m qubits, simulation of 
one fermionic gate takes 0(m) qubit gates and vice versa. We show that using different 
encodings, the simulation cost can be reduced to O(logm) and a constant, respectively. 
. Nearest-neighbors fermionic gates on a graph of bounded degree can be simulated at a 

constant cost. A universal set of fermionic gates is found. We also study computation 
with Majorana fermions which are basically halves of LFMs. Some connection to qubit 
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quantum codes is made. 



Introduction 

The notion of locality plays the key role in the definition of computation process. The same basic 
principles apply to classical computers and the circuit model of quantum computation [jXJ, 

1. The computer consists of small pieces, or cells (bits, qubits, qutrits or something else). 

2. It is allowed to operate on few cells at a time. 

3. All cells are identical, so each operation has a model which can be applied to different 
sets of cells. We call such a model a gate (like the CNOT gate) while the operation itself 
is called a gate application (like CNOT applied to qubits 5 and 8). 

*On leave from L. D. Landau Institute for Theoretical Physics 
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In fact, the main difference between classical computation and quantum computation is the 
concrete interpretation of these postulates. The standard quantum interpretation is as follows. 

li. Each cell is described by a Hilbert space of small dimensionality. (Without loss of gen- 
erality, these spaces are two-dimensional, in which case the cells are called qubits). The 
Hilbert space of the entire computer is the tensor product of the spaces associated to the 
individual cells. 

2i. Each operation is described by a unitary operator which is the tensor product of some 
operator U, acting on the selected qubits, and the identity operator acting on the rest of 
the system. 

3i. A p-qubit gate can be defined as U acting on some standard 2 p -dimensional space (which 
is the same as a standard set of p qubits). 

This interpretation might be perfect from the complexity-theoretic point of view, but it 
does not necessarily correspond to physics. At the fundamental level, Fermi systems do not 
satisfy the condition (2i). Hence using fermions as carriers of quantum information H should 
be considered as a different computation model, although it is equivalent to the standard one in 
a certain sense. At the macroscopic level, collective quantum degrees of freedom (or excitations, 
such as anyons 0, |[) do not even satisfy the condition (li). In all such cases, it makes sense to 
abstract the nontrivial locality properties from physical details. This will lead us to a definition 
of "fermionic gates", "anyonic gates" and other quantum computation models which deserve 
careful study. 

What can one expect from alternative models of quantum computation? It is very unlikely 
that any physical system would provide more computational power than the standard quantum 
computation model has. (This might be wrong for quantum gravity but here we can only guess). 
So, the alternative models should be polynomially equivalent to the standard one. There are 
indeed several results which support this statement. Firstly, D. Abrams and S. Lloyd || have 
shown that a system of m local fermionic modes (i.e. sites which can be empty or occupied) 
can be simulated on a quantum computer in such a way that one fermionic gate takes 0(m) 
qubit operations. In some cases (one of the assumptions is that the number of particles is 
conserved) faster simulation is possible. We will extend this result by showing that in the 
general case (when the number of fermions is conserved only modulo 2) each fermionic gate can 
be simulated by O(logm) qubit gates. Secondly, TQFT computation (which is more general 
than quantum computation with anyons) can be simulated in polynomial time on an ordinary 
quantum computer 0. 

Thus alternative quantum computation models do not generate new computational classes. 
Rather, they provide new descriptions of the standard class BQP (the class of problems that 
are solvable on a quantum computer in polynomial time). These descriptions may be useful 
to find new quantum algorithms, error-correcting codes, fault-tolerant procedures, or to prove 
that BQP is contained in some other computational classes. They may also open the door to 
new physical implementations of a quantum computer. 

1 A more general notion of locality 

Why isn't the standard interpretation of locality good in all cases? The answer is in the way 
we describe quantum evolution. The Hilbert space and state vectors are very convenient tools 



2 



but they are not directly related to physically observable things. Operators are much more 
"real" as they allow one to describe interaction between the system and the rest of the world. 
In fact, the causality principle is usually stated in terms of operators: "operators in spatially 
separated points commute". Our definition of locality will be in the same spirit. 

A quantum system can be characterized by a finite-dimensional C*-algebra[] Q whose ele- 
ments are called "physical operators". As a matter of fact, they are operators on a suitable 
Hilbert space. Indeed, it is a theorem that Q can be represented as 0jL(W 3 -), where L(£) 
stands for the algebra of operators on the space C. Hence Q acts on Tt = QjTCj. (In the case 
of fermions, H = Ho © H\ is the Foch space split into the subspaces corresponding to an even 
and odd number of particles; Q is the algebra of operators which preserve the parity). 

To define locality, we will assume that the system is associated with some set of sites M. 
The following properties are postulated: 

1. For each subset of sites S C M, there is a C*-subalgebra Q(S) C Q. Elements of Q{S) are 
called physical operators acting on S. We require that Q{M) = Q, C/(0) = C • I (where 
C is the algebra of complex numbers, / is the unit element of Q), and G(S) C Q(S') if 
S C S'. 

2. If Si H S2 = then any two operators X\ G G(S\) and X 2 G commute. 

The concept of unitarity is well defined in this setting: an element U G Q{M) is called 
unitary if UW = WU = I (the operation f is a part of the C*-algebra structure). Note that 
nonunitary elements of Q(M) also have physical meaning because they can be used to construct 
a unitary operator on a larger space H® C, where £ represents some external system (e.g. a 
measurement device). Such an operator generally has the form U = Y^k Ak® Bk (A^ G Q(M)), 
i.e. U G Q{M) ®L(C). 

Thus we have given a more general interpretation of the locality postulates 1 and 2 which 
were discussed in the introduction. (We put aside the postulate 3 here). 

2 Local fermionic modes 

Consider m sites (numbered through m— 1) each of which can be either empty or occupied by a 
spinless fermionic particle. Such sites will be called local fermionic modes (LFMs). The Hilbert 
space H of this system, known as Foch space, is spanned by 2 m basis vectors |n , . . . ,n m _i), 
where rij = 0, 1 is the occupation number of the j-th site. Everything related to fermions can 
be expressed in terms of annihilation and creation operators a-poj, (j = 0, . . . ,m — 1). The 
operator a,j acts on basis vectors as follows: 

a j \n ,...,n j „i,l,n j+1 ,...,n m - 1 ) = (-l)^==o ns |n , . . . , n 3 -_ l3 0, n j+1 , . . . , n m _i), ^ 
a j \no,...,n j -i,0,n j+1 ,...,n m - 1 ) = 0; 

a,j is the Hermitian conjugate. Note that the definition depends on the order of LFMs! (One 
may rather say that the basis depends on the order while aj,aj do not, since all relations 

1 C*-algebra is a generalization of the algebra of linear operators on a Hilbert space. The properties of 
operator addition, multiplication, Hermitian conjugation and the operator norm are axiomatized in a certain 
way. However, instead of using the axioms, we will rely on a characterization of finite-dimensional C*-algebras 
given below. 
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between them are permutation- invariant). The annihilation and creation operators generate 
the algebra T = L(7i) and have the following commutation rules: 

0, 

0, (2) 
Sjk- 

The Hilbert space of m LFMs splits into two parts: Ji — Tio®TCi, where "0" and "1" refers 
to the total fermionic parity n = J2Y=o n i ( m °d 2). Physical operators are those which preserve 
the parity. Note that the Hamiltonian of a real Fermi system always satisfies this condition] 2 ] 
unlike the operators %,aj alone. The algebra of physical operators T = L(7^o) © L(Wi) is 
spanned by products of even number of aj, a]-. (The notation Q in the previous section referred 
to the general case whereas T is reserved for LFMs). 

Let SC. {0, . . . , m — 1} be a set of LFMs. Physical operators on S are linear combinations 
of even products of aj,Oj, j G S. The algebra of such operators is J-'iS) = F(S) fl J 7 , where 
F(S) C T is generated by aj, al, j G S. The conditions (1) and (2) from Sec. [TJ are obviously 
satisfied. Moreover, J-"(Si) commutes with F(S 2 ) if SiH S 2 = 0. 

The occupation number rij is an eigenvalue of the operator hj = ajeij G ^F({j}), which 
means it can be measured locally (by acting on the j-th LFM and some external device). The 
occupation number can not be changed locally, though. 

Here are some examples of unitary operators acting on one or two LFMs: exp^/fajoj) (action 
by an external potential), exp(ifla\aja\.ak) (two-particle's interaction), exp(z(7ajafc + j*a\aj)) 
(tunneling) and exp(i(jakaj +7*aja£)) (interaction with a superconductor). We will show that 
these operators (for all or for some particular values of (3 G R and 7 G C) form a universal 
set, i.e. any unitary operator can be represented as a composition of these ones to any given 
precision, using ancillas. 

In our computation model we allow to use ancillas in the state |0). (To be more accurate, 
we should say that the input state is padded by some number of zeros to the right. Actually, 
the order does not matter in this case). If we speak about implementing a unitary operator, the 
ancillas must return to the state |0) by the end of the procedure. As is usual, this restriction 
does not apply to computing a Boolean function which proceeds as follows. One starts from 
a basis vector \n , . . . , n m _i) representing the input data, adds some ancillas, applies some 
sequence of local unitary operators and reads the result by measuring some of the occupation 
numbers. 

3 Relation between LFMs and qubits 

The Hilbert space of m LFMs can be identified with the Hilbert space of m qubits B® m (where 
B stands for the two-dimensional space C 2 endowed with the standard basis {|0), |1)} ): 

|n ,ni, . . . ,ra m -i) = |n ) ® \ni) <g> • • • |n TO _i), 77^ = 0,1. (3) 

2 In electon systems, the Hamiltonian also preserves the electric charge, so terms with different numbers of 
aj and aj are usually forbidden. Our model is mostly relevant to superconductors where the total charge of 
excitations is not conserved, so terms like ajOk appear in the effective Hamiltonian. 
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Measurement of rij is the same as eigenvalue measurement of cr|. A physical fermionic oper- 
ator corresponds to a qubit operator which preserves the parity, i.e. commutes with Ilj^o 1 a ]- 
However, "applying a gate to a set of LFMs" is very different from "applying a gate to a set of 
qubits" . 

Let X be a parity-preserving p-qubit operator acting on qubits numbered through p — 1. 
Applying it to the qubits j , . . . ,j p -i is a straightforward procedure. The Hilbert space of m 
qubits E® m can be identified with B m <g> £?®( m -p) by the qubit permutation P : \n , . . . n m _i) h-> 
|n, , . . . n J - p _ 1 ) ® | the other n,-). Then the action of X is defined as follows: 

X[j ,...,jp-i] = P~\X ® I m(m - p) )P. 

If we want to apply X to the LFMs jo, ■ ■ ■ ,jp-i, the procedure is different. First, we 
should expand X into products of ao, . . . , a p -i, aj, . . . , aj^. Then we replace each a r by a Jr 
and each a\ by aj r . The resulting operator will be denoted by X{j , . . . ,j p _ 1 }. For example, if 
X — |1, 0)(0, 1| = ajai then Xjj, A;} = ajja fc . This operator acts as follows: 

a]a k \...,0,n j+1 ,...,n k - 1 ,0,...) = 0, 

okfc | . . . ,0,n j+ i, . . . ,n fc _i, 1, . . .) = (-1)^.=j+i ™ 3 | . . ., l,n i+ i, . . . ,n fc _i,0, . . .), 

a)a k | • • • , l,n j+1 , . . . ,n fc _i,0, . . .) = 0, 

a)a k | l,n j+1 , ... 1, .. .) = 0. 

Not only X{j, A;} 7^ X[j, fc] but also fc} is non-local in terms of qubits: it involves all the 
qubits with numbers from j to k. 

A unitary qubit gate (LFM, or fermionic gate) is a unitary operator U meant to be applied 
to a number of qubits (LFMs); a p-qubit or a p-LFM gate acts on the standard Hilbert space 
B® p . Operators of the form U[j,k] or U{j,k} are called gate applications. We will usually 
consider unitary gates up to overall phase factors. A set of gates is also called a basis. A circuit 
of size s in a basis A is a composition of s applications of gates from A, i.e. an expression 
of the form U s {j Sy0 , . . .,j s , Ps -i} ■ ■ ■ ^i{ji,o, • • - ,ji, P1 -i}, where U k G A. Such an expression is 
evaluated by a unitary operator which is said to be represented by the circuit. 

Note that X{j} = X[j], so one-LFM gates are simply parity-preserving one-qubit gates. 
Up to an overall phase, such gates have the form A(e^), where A(U) denotes the controlled U . 
(If U acts on p qubits then A(U) acts onp + 1 qubits; in our case p — 0). 

More generally, X{j, . . . , j+p— 1} = X\j, . . . , j+p— 1]. This allows one to represent 
fermionic gates in terms of qubit gates and vice versa. We will now show how to do that in the 
case p — 2. 

Suppose we want to execute a two-LFM operator X{j, k} (w. 1. .0. g. j < k). First, we move 
the A:-th qubit next to the j-th one by swapping it with its nearest neighbors. Then we apply 
X{j,j + 1} = X[j,j + 1\ and move the A;-th qubit back to its original position. However, what 
we actually need here is to interchange LFMs, not qubits. This is different even if the LFMs 
(qubits) are next to each other! 

A swap between two qubits (with numbers and 1) is defined in the obvious way: (<->) : 
\ n o, n i) l— *■ l^ij^o)- A swap between two LFMs is a unitary operator (•<=>■) such that 

(^)a (^) f = ai, a 1 (^) f = a . (4) 
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These equations have a unique solution (up to an overall phase factor): 

(■&) = I - ala - a\ai + a\a + ajai = exjp(i^-(al - a\)(a - aij) : 
This differs from the qubit swap by the "-" sign. So, 



0,0) - 


* |o, 0> 


0,1)- 


*|i,o) 


1,0)- 


Ho,i) 


1,1)- 


■> — 11, 1 



(5) 



(6) 



where D = A(a z ) : |a, 6) i— > (— l) a6 |a, 6) is a "swap defect" operator. 

To perform the procedure described above, we do not have to actually interchange LFMs 
or qubits. Instead, we can simply apply operators D{l,r} = D[l,r] to all pairs of qubits 
that otherwise would be interchanged. So, any two-LFM operator X{j, k} (j < k) can be 
represented by a qubit circuit as follows: 

X{j,k} = D[k-l,k]---D\j+l,k]X\j,k}D\j + l,k}---D[k-l,k}. (7) 

Conversely, any parity- preserving two-qubit operator is represented by a fermionic circuit: 

X\j,k] = D{k-l,k}---D{j+l,k}X{j,k}D{j + l,k}---D{k-l,k}. (8) 

This method also works for operators which act on more than two LFMs (qubits). 



4 A universal set of LFM gates. 

We have shown that fermionic gates are equivalent to parity-preserving qubit gates modulo the 
swap defect operator D. So, to obtain a universal set of fermionic gates, we only need to find 
a universal set of parity-preserving qubit gates. It will be possible to represent the operator D 
by these gates too, exactly or approximately. (If only approximate representation is possible, 
one may have to pay extra cost when simulating a single fermionic gate because the more qubit 
gates are used, the more accurate they should be. We will avoid this problem, though). 

We claim that the following gates are sufficient to represent any parity-preserving operator 
to any given precision (using ancillas): 

A(e^ 4 ), A(a z ), H : \a, b) - -L £(_i)*| a +&+ C| c). (9) 

(Here a, b, c E F 2 = {0, 1}, so the expression a + b + c is considered modulo 2). Note that 
D = A(a z ) belongs to this set of gates. 

Any parity-preserving operator U can be considered as a pair of operators (Uq, Lq), where 
U acts on the even sector Ho whereas U\ acts on the odd sector Hi. We will first show how 
to get a given U while not caring about Lq — we will actually get Lq = U - (The comparison 
between U and Lq is made through identifying Ho and H\ by the map (^[O]). Then we will 
implement operators of the form (/, Y) which can be used to correct the first step. 

Any operator X on m — 1 qubits can be transformed to a parity-preserving operator X on 
m qubits by using the extra qubit to maintain the parity: 

X = V' 1 (I B ® X) V, V : \n Q ,ni, . . . ,n m -i) -> \n Q +. . .+ra m _i, ni, . . . ,n m -i). (10) 
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Note that the gate H (see eq. (|9|)) corresponds to the Hadamard gate H = ^= y ^ ^ J . If 

X already preserves the parity, as it is the case with the operators A(e i7r//4 ) and A(<r z ), then 
X = Iq®X. (This equality is actually a characteristic property of parity-preserving operators). 

The operator V is unitary. It maps the even sector TCq onto the subspace C which consists 
of vectors |0) ® |£) G Any operator on this subspace extends to an operator of 

the form ig <g> X. Hence any Uq G L(Wo) extends to an operator of the form X. Note that 
V and Ib® X commute with a x [0], so X also commutes with a x [0}. Therefore X = (Uq,Ui), 
where U\ = a x [0] Uqct x [0}, or simply U\ = Uq if the identification between TCq and Hi is used. 

If an operator X is represented by a quantum circuit Al ■ ■ ■ Ai, one can replace each gate 
application A^ by A^ to get a quantum circuit for X. Indeed, eq. fllCf ) defines a *-algebra 
homomorphism, i.e. 

x^+x 2 = T 1 +T 2 , ^x = cx, x^x 2 = x 1 x 2 , 16 = x\ T=i. 

It follows that any universal gate set A transforms to a set of parity-preserving gates A which 
are universal on the even sector. The following gate set is known to be universal [H: A = 
|A(e l7r / 4 ), A(a z ), H\. The corresponding gate set A is given by eq. (0). (The parity-preserving 

gates A(e l?r//4 ) and A(a z ) are copied from ^4 to A unchanged). Thus the basis @ allows one to 
obtain at least unitary gates of the form (Uq, Uq). 

With parity-preserving operators, we can still use one of the standard techniques in quantum 
circuit design — gates with quantum control. Indeed, if U preserves the parity then A(U) also 
does. More generally, A(X) = A(X) (up to a permutation of the control qubit and the parity 
qubit). So, if we can represent A(X) by a circuit,]] we can also obtain a circuit for A(X) by 
the procedure described above. For example, the operator A(a 1 ')[l,2] can be represented as 
H[2] A«)[l,2]if[2], hence 

A(^)[l,0,2] = H[0,2] A(a 2 )[l,2] #[0,2]. 

Here we use qubit to maintain the parity whereas 1 is the control qubit. (The notation 
A(X)[. . .] suggests that the control qubit goes first). 

Now we are in a position to implement operators of the form (J, Y) using the gates @. We 
will also use one ancilla which will be assigned the number m. First, we execute the operator 

W: \rio,n 1 ,...,n m - 1 ,n m ) h-> \n + . . .+n m - 1 , m, . . . , n m _i, m + . . .+n m ), 

W = A(^)[m-l,0,m] ■•• A(^)[l, 0, m], ' 

where A(a x ) : \a,b,c) i— > \a, b+a, c+a). Now qubit indicates the total parity. 
Let X = (Y, Y). Since X can be represented according to (|10D , and because 

W = V^m^, . . . ,m-l] V, where V = V[0, . . . ,m—l), 

3 Implementing A(X) in the basis A is only slightly harder than implementing X. Indeed, for each gate X 
from the basis A, the operator A(X) can be represented by a circuit in the same basis exactly. Therefore, the 
circuit for A(X) will be larger than the circuit for X only by a constant factor. Note, however, that X may 
have been implemented up to a phase factor; this phase factor becomes important when we consider A(X). The 
necessary correction can be achieved by an operator A(e I< ^) which should be implemented separately. 
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we conclude that 



W- 1 A{X)[0,m,l,...,m-l]W = V' 1 A(X)[0, 1, . . . , m-l] V = U®I B , U=(I,Y). (12) 



Note that in this case, the ancilla is not affected no matter what its initial state was; one can 
even use a data qubit as the ancilla. 

We have proved that the gate set @ is universal in the class of parity-preserving operators. 
It remains to represent these gates in terms of creation and annihilation operators. The first 
two are simple: 

A(e l7r / 4 ) = exp(z— ajao) , A(a z ) = exp(i7ralaoa\ai). (13) 



Unfortunately, the gate H in the fermionic representation looks ugly, so we first represent it as 
follows: 

#[0,1] = A(-z)[l]G[0,l]AH)[l], ' 1 1 ' 



G 



i 1 



So, the fermionic gates (0), together with 



G = exp(— i-^(a — <4)(ai + a\fj = exp(i—(alai + a[a )^j exp(i—(aia + ala\) ) > (■'■■!) 
form a universal set. Obviously, this gate set is also universal: 



exp^i|(aQai + a\a 



exp(ij(aia + aja|)y. 
exp(z7raoaoalai) 



> . 



(15) 



5 Fast simulation procedures 

So far we have been using the standard identification (J3|) between the Foch space 7i and the 
Hilbert space of m qubits B® m '. This identification has allowed us to consider qubit gate 
applications X[j,k] e L(i3 (X " n ) and LFM gate applications X{j, k} G T C L(TC) as operators 
acting on the same space. Now we are going to discuss a more general way of simulating LFMs 
by qubits (or vice versa). Let J be an encoding of m LFMs by m! qubits, i.e. J : H — » B® m 
is a unitary embedding. (J being a unitary embedding means that J^J = In- Note that J,P 
is the projector onto C = Im J C £>® m ). We say that an operator U' G L(i3® m ) represents an 
operator [/ G L(W) if 
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| J commutes, 



i.e. JC7 = U'J. 



(16) 



The operator U, as well as the LFMs or qubits equivalent to them, are called logical. The mf 
qubits the operator U' acts on are called code qubits. 

In this section we will show that each LFM gate can be simulated by O(logm) qubit gates. 
Surprisingly, this does not require quantum codes in the proper sense, i.e. J : Tt — > B® m is a 



8 



map onto. Moreover, J takes basis vectors to basis vectors. We will also show how to simulate 
qubit gates by LFM gates at a constant cost using a subspace of the Foch space. 

Gate action on a set of qubits (LFMs) can be described as follows. First, we extract two 
(or some other number) qubits or LFMs from the quantum memory which is now considered 
as a "black box". We place these qubits at positions —2 and —1. Then we apply the gate and 
put the qubits (LFMs) back into the memory. In this description, it does not matter what gate 
we apply, a qubit one or a fermionic one (because X{ — 2, —1} = X[— 2, — 1] ). Once extracted, 
LFMs can be regarded as qubits — the difference is in the way we extract them. With qubits, 
we apply the operator 

[<-j] : B® m — > B®B® m : \n^)®---®\n d )®---®\n m _ x ) ^ \n 3 ) ® |n )<8)- • -®|0)(8)- • -<g> 

This is a unitary embedding into a larger Hilbert space; it can be represented as adding an 
ancilla in the state |0) followed by a unitary operator. With fermions, we move the j-th LFM 
through all the LFMs left to it, so we should take swap defects into account. Thus we get 
another unitary embedding: 

: n^B®n, 

{^j} \ n o, . ,n m _i) = (-lp^ s =o ns \ nj ) g) | no? . . . ,0, . . . ,n m _i). 

It is easy to verify that our recipe is correct, i.e. applying the operator X[— 2, —1] after extracting 
two LFMs is equivalent to applying X{j, k} before the extraction: 

X[-2,-l](l B ®{^k}){^j} = (l B ®{^k}){^j}X{j,k}. (18) 

(Our notations are somewhat confusing so this comment could be helpful. When we extract 
the j-th LFM by applying {^=j}, we add one qubit at position —1. When we then extract the 
k-th LFM, this qubit moves to position —2 while another one is being inserted; this can be 
described by the operator I& ® {<=k}.) 

Suppose we want to simulate LFMs by qubits at low cost. The problem with the standard 

encoding (0) is that multiplication by the factor (— l) nj ^s=o Us requires too many operations. 
The simplest solution would be to store yj = Y? a =o n s instead of rij. (Remember that we consider 
n s as residues (mod 2), so the sum is also taken (mod 2) ). This does solve the problem but 
also creates a new one: when rij becomes as a result of extraction, we have to modify all yk, 
k > j. So, we need to balance the complexity of computing J2lZo n s and that of updating the 
encoded quantum memory when rij changes. This kind of trade-off can be achieved by storing 
some partial sums Y J h s = a n s- 

In general, we will use encodings of the form 

Xj = n si S U) C {0, . . . ,m - 1}. 
We start with an example of such an encoding for m = 8 (the diagram next to the equation 
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illustrates grouping of n s into Xj 

n 2 x 



n ni 



n 2 n 3 



n 4 n 5 



n 6 n 7 



xq = n X2 = n% x 4 = 714 %e> — n§ 
x i = n o + Tii x 5 = ri4 + n 5 

x 3 = n + ri! + n 2 + n 3 
x 7 = n + rii + n 2 + n 3 + 714 + n 5 + n 6 + n 7 

A binary tree structure is apparent here. To proceed, we will represent the qubit indices 
0, . . . ,m — 1 by binary strings. The length of these strings is not fi xed; w e may add an arbitrary 
number of zeros to the beginning of a string, e.g. 3 = 11 = Oil = 0011. 

Let us define a partial order ^ on the set of binary strings. We write a t -i . . . «o ^ Pt-i • ■ ■ A) 
if cxi = A f° r I > while A -i = . . . = (3q = 1 (for some Z )- F° r example, 



000 -< 001 
010 



-< 011 



> -< 111, 



100 -< 101 



110 



where j -< k means that j ^ k but j ^ k. Note that if j -< k then j < k. 
Now we can specify our encoding for arbitrary m: 



\n , . . . ,n m _i) i-> |x ) <g> • • • ® |x m _i), 



E 



».s 



(19) 



It is important that each n s enters only O(logm) of Xj. Hence applying a x [s] to one logical 
qubit amounts to logarithmically many a x gates being applied to the code qubits. 
The inverse transformation (from Xj to n s ) is also simple: 



x. 



Xs i 



a t -i . . . «o G K(3 t -i ■ ■ ■ A>) if and only if 

Pi = . . . = Po = 1, ati = 0i for / ^ l , ai 

The sum in ( f20|) contains only O(logm) terms. Moreover, yj 
as a sum of O(logm) numbers x s : 

Uj = E Xs > 

a t -\ . . . ao E L^A-i . . . ~/3qJ if and only if 

A = 1, on = A for / > Zq, a io = 0, a io _i = . . . = 



(20) 

(for some Iq). 

Y^l=o n s can be also expressed 



(21) 



a = 1 (for some Z 



It remains to actually represent the LFM extraction operator {<t=j} by acting on the code 
qubits. Since {^j} increases the number of qubits by 1, we should add this extra qubit first. 
We place it at position —1 and initialize by |0). Then we make its value equal to rij by applying 
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the operators A(a x )[s, — 1], s G K(j) U {j} (see eq. (pD[)). After that, can be turned into 
by executing the operators A(a x )[— 1, fc], fc ^ j (see eq. (0)). Finally, we multiply by (— \) n i y i 
by applying A(a z )[— 1, s], s G L(j) (see eq. (|2l|) ). To summarize, we execute the operator 

U'= nA(0[-MlIlAK)[-U] II AM[«-1] (22) 

which has the property (ijg ® J) {^j} |£) = £/'(j0) £g> J|£)) f° r an y 10 e This requires 
O(logm) operations. 

Simulating qubits by LFMs is easier and faster. One can use this simple encoding: 

|n ) ® |ni) <g> • • • ® |n m _!) h-> |n ,n 0) ni,m, . . . ,n m _i,n m _i), (23) 

i.e. each qubit is represented by a pair of LFMs with an even number of fermions, |0) i— > 1 00) , 
|1) i — ► 1 11). Each two-qubit operator X is represented by a four-qubit operator X' . The 
operators X'[2j, 2j + l, 2k, 2k + 1] and X'{2j,2j + l,2k,2k + l} act the same way on the code 
subspace. So, when we use the encoding (^3|), it does not matter whether the code elements 
are qubits or LFMs. The simulation cost is just a constant. 



6 Majorana fermions 

It is possible, at least mathematically, to split each local fermionic mode into two objects. 
These halves of LFMs are called Majorana fermions.0 

Let us introduce a set of Hermitian operators Cj (j = 0, . . . , 2m — 1): 

k—l t k— 1 

c 2k = a k + a{ = a x [k] JJ o*\j\ , c 2k+1 = = a y [k] ^[j]. (24) 

j=0 1 j=0 

These operators satisfy the commutation relations 

CjC k + c k Cj = 25 jk . (25) 

We can define new locality rules on the algebra T = L(7Y ) © L(^i) which will be now denoted 
by M.. (Also T = L(7i) will be denoted by M). We say that there are 2m sites called Majorana 
fermions. For each set of sites S C {0, . . . , 2m — 1}, let M.(S) C M. be the subalgebra generated 
by Cj, j G S\ (Such an algebra is known as complex Clifford algebra). Then M.(S) = Ai(S)nA4, 
i. e. physical operators on S are linear combinations of even products of Cj, j G S. 

According to this definition, Majorana fermions "exist" in any Fermi system. A nontrivial 
thing is that it is possible (at least theoretically) to pair them up by interaction, so that few 
Majorana fermions remain unpaired and separated from each other 0. Such systems could be 
used as decoherence-free quantum memory. Indeed, a single Majorana fermion can not interact 
with the environment by itself (because the operator Cj is not physical), so the decoherence 
can only arise from environment-mediated interaction of two Majorana fermions. But if they 
are well separated in space, such interaction should be exponentially small (a finite correlation 
length in the environment is assumed). Roughly speaking, we keep two halves of a qubit apart, 
so the qubit is decoherence-free! 

4 In field theory, this term usually means something more specific, but it is sometimes used in this sense too. 
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An application of a Majorana gate is denned as in the case of LFMs. Let X be a physical 
operator acting on p Majorana fermions numbered through p — 1. Then X can be expanded 
into (even) products of Cj, j — 0, . . . ,p— 1. If we substitute Cj r for c r , we will get an operator 
A{{j , • • • ,j P -i}} e M({j , . . -Jp-i})- 



To find a universal set of Majorana gates, it suffice to rewrite the operators (|13[), (|14|) in 
terms of cq = cio + aj, ci = — i(ao — aj), c 2 = ai + a\, c 3 = — i(a\ — a\). 

aja = |(1 + icoci), a|ai = |(1 + zc 2 c 3 ), 
exp(^aja ) = e m/8 exp(-|c ci) , exp(-z|(a _ + = exp(fcic 2 ), 

exp(z7raja alai) = e" m/4 exp(f c Ci) exp(fc 2 c 3 ) exp(if c cic 2 c 3 ) . 
Hence the following gate set is universal (up to phase factors) : 

jexp(^coci), exp(^c cic 2 c 3 ) j . (26) 



The second gate in the set ( p6|) describes four-particle interaction. When it comes to physical 
implementation, this gate will be particularly difficult to realize. Unfortunately, it is not possible 
to do universal quantum computation by acting on three or fewer Majorana fermions at a time. 
This way one can only generate the group of unitary operators which act by conjugation as 
follows: 

U Cj rf = J2PjkC k , (27) 
k 

where (3 G SO(m), i.e. (fljk) is a real orthogonal matrix with determinant +1. 

7 An alternative to the four-particle Majorana gate 

In this section we show that, for the purpose of universal computation with Majorana fermions, 
the gate exp(z|c cic 2 c 3 ) can be replaced by a nondestructive eigenvalue measurement of the 
operator CoCiC 2 c 3 . (A measurement being nondestructive means that vectors in each of the 
eigenspaces remain intact; in other words, no extra information is learned or leaks to the 
environment). Such measurements might be easier to implement, they are also useful in some 
theoretical application of fermionic computation [|1(| . 



Let us assume that the following operations are possible: 

1. Applying the unitary gate R = exp(|coCi). (Note that RcqB) = — ci, Rc\R} = cq). 

2. Creation of an ancilla pair in a state which is the eigenstate of c 2 fcC 2 fc + i = i(l — 2a\at) 
corresponding to the eigenvalue i. 

3. Eigenvalue measurement of CjCk- 

4. Nondestructive eigenvalue measurement of CjCkC r c s . 
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Moreover, we can base the choice of the next operation on the previous measurement outcomes. 
(One may call such quantum computation adaptive). Of course, the amount of classical com- 
putation involved in this choice should not be too large, so we should better include it into the 
overall size of the quantum circuit. 

Remark. In the standard scheme of quantum computation, measurements in the middle of 
computation are redundant and can be simulated by unitary gates (e.g. see |TT[). However, 
this is only true if one uses a universal set of unitary gates. Otherwise measurements and 
adaptiveness can indeed add extra power to unitary operators. 

Suppose we want to apply the operator exp(i|coCiC 2 c 3 ). Let Majorana fermions 4 and 5 
form an ancilla pair, so the input state of the system satisfies 

(c 4 + z C5 )|^ in ) = 0. (28) 

We measure the eigenvalue of C0C1C3C4. The outcome is either +1 or —1, which means that |^i n ) 
gets multiplied by the projector U+\ = ~(1 + c cic 3 c 4 ) or II if = |(1 — c cic 3 c 4 ), respectively. 

(More exactly, \^- m ) >— ► P±i^ 2 n±]|\I/hi), where p z is the probability to get outcome z). Then we 
measure the eigenvalue of C2C4. The possible eigenvalues +i and — i correspond to the projectors 
Il+i = |(1 — zc 2 c 4 ) and H^j = ~(1 + ic 2 c 4 ). We claim that after some correction depending 
on the measurements outcomes, we will effectively execute the operator exp(i|coCiC2C 3 ) while 
leaving the ancilla pair intact. Indeed, 

exp(if C0C1C2C3) |^ in ) = 

= 2 exp(f c 2 c 5 ) i(l - ic 2 c 4 ) |(l + c cic 3 c 4 ) \% n ) = 
= 2iexp(fc ci) exp(fc 2 c 3 ) exp(f c 2 c 5 ) \{l - ic 2 c 4 ) |(l - c cic 3 c 4 ) |* in ) = (29) 
= 2iexp(f c ci) exp(f c 2 c 3 ) exp(-~ c 2 c 5 ) |(l + ic 2 c 4 ) |(l + c cic 3 c 4 ) |* in ) = 

= 2 exp(-f c 2 c 5 ) |(l + ic 2 Cij \ {l - c cic 3 c 4 ) |^ in ) 

(we have used eq. (p8|)). Thus we can apply a suitable correction operator U yz in each of the 
four cases (C/ +i)+ i = exp(|c 2 c 5 ) if the outcomes were +1 and +i, etc.) so that 

exp^cocipaca) |* in ) = 2 U yz U^ \*- m ). 

Each of the four outcome combinations occurs with probability 2 -2 = j. The final state is 
always the desired one, \^/sn) = exp(ijCoC\C2C^j \^- m ). 

8 Superfast simulation of fermions on a graph 

The results of Sec. [| suggest that fermions have slightly more computational power than qubits. 
The logarithmic slowdown in simulation of fermions seems to be inevitable in the general case. 
However, in the physical world fermions (e. g. electrons) interact locally not only in the sense 
that the interaction is pairwise, but also in the geometric sense: a particle can not instantly 
jump to another position far away. Such physical interactions might be easier to simulate. In 
this section we study an abstract model of geometrically local interactions. The result is that 



13 



geometrically local gates can indeed be simulated without any substantial slowdown, i. e. the 
simulation cost is constant. Therefore one can speculate that, in principle, electrons might not 
be fundamental particles but, rather, excitations in a (nonperturbative) system bosons. Of 
course, this is only a logical possibility which may or may not be true. 

What follows is a definition of the model. Consider a connected unoriented graph F = 
(M, E), where M — {0, . . . , m — 1} is the set of vertices, and E C M x M is the set of edges. 
(As the graph is unoriented, (j, k) and (k, j) either both belong to E or both do not belong to 
E) . We will assume that the degree of each vertex is bounded by some constant d. Let us put 
a local fermionic mode on each vertex. We will consider only the even sector of the system, 7i , 
i. e. the total number of fermions is required to be even. The allowed unitary operations are 
one-LFM gates and two-LFM gates applied to any pair of vertices connected by an edge. 

We are going to identify the Hilbert space 7i with a codespace of a certain symplectic 
(stabilizer) code so that each elementary operator of the form a\a k , as well as aja fc , a\a^ 

ajak or a k aj, where (j, k) G E, be represented by operators acting on 0(d) qubits. Then each 
one-LFM gate and each two-LFM gate applied to neighboring vertices will be also represented 
by an operator acting on 0(d) qubits. As d = const, this means one can simulate each of the 
allowed fermionic operations by a constant number of one-qubit and two-qubit gates. 

It will be convenient to use the Majorana fermions operators C2k, C2k+i (see eq. ([24])) instead 
of a k , a\. The list of elementary operators to be represented can be reduced to these ones: 



B k = —ic 2 kC2k+i for each vertex k, 
Ajk = —ic2jC2k for each edge (j, k) E E. 



(30) 



These operators satisfy the following relations: 

Bl = B k , A] k = A jk , B 2 k = l, A% = 1, A kj = -A jk , (31) 

B k B t = B t B k , A jk B t = (-l) 5 ^ 5kl B l A jk , A jk A ls = (-l) 5 ^ +s ^ A ls A jk , (32) 

i p 'Ajo^Afrjz ■ ■ ■ Aj p _ 2 j x Aj u j = 1 for any closed path on the graph. (33) 

It is easy to prove that B k , Aj k modulo these relations generate the algebra of physical operators 
T = L(7io) © L(7Yi). However, we are considering only the even sector now. Having been 
restricted to Tio, the operators B k satisfy an additional relation (which was false in T): 

\[B k = l. (34) 
k 



Hence the algebra L(7i ) is generated by B k , Aj k modulo the relations (pl|)- (|34]) . 

To construct the code, we put a qubit on each edge of the graph. Thus a^ k = o~ k j denotes 
the Pauli operator a a (a = x,y,z) acting on the edge (j,k). The operators B k , Aj k defined 
above will be identified with some operators B k , Aj k acting on the code subspace C (which will 
be defined later). We start with defining the action of B k and Aj k on the entire Hilbert space of 
the qubits. Our construction depends on two arbitrary choices. Firstly, we choose orientation 
for each edge of the graph. This can be described by a matrix (e^) such that e k j = —€j k , 

5 The term "stabilizer code" has become traditional but it is somewhat confusing because any code can 
be defined in terms of stabilizer operators. (It is actually practical to do so for nonbinary codes related to 
anyons §). We prefer to use the more specific terms "symplectic code" . 
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€jk = ±1 for each edge (J, k) G E. Secondly, for each vertex k, we order all incident edges 
(j, k). This order will be denoted by <. Now we put 



k 



B k = n °i» 

3- (j,k)£E 



Ajk = tjkVjk n °ij n 



°sk- 



(35) 



l.(l,j)<(k,j) s:(s,k)<(j,k) 

j k 



These operators satisfy the relations analogous to (|31~1) , (j32|) and (|4])> but not (p3| ). 

Finally, we define the code subspace C C B® u (where u is the number of qubits) by imposing 
stabilizer conditions: 

|V>) G £ if and only if (5,- Or ..,i JJ _ 1 |'0) = |t/>) for any closed path (j , . . .Jp-i,jo), 

(36) 

a * . - apa. . A A ■ ■ A- 

^30,-,3p-1 L ^30,31^31,32 ^p-Uv-l 3p-l,30 m 

The stabilizer operators C , J - 0i ...j J) _ 1 are Hermitian and can be represented in the form ± riy.fe) ^a^j 
where (j, k) runs over a set of different qubits. The set of stabilizer operators is obviously 
redundant but it is consistent, meaning that (i) they commute with each other, and (ii) whenever 
the product of several stabilizer operators is a constant, this constant is 1. The number of 
independent stabilizer operators equals the number of linearly (mod 2) independent cycles 
which in turn equals u — m + 1. Hence 

dim £ = 2"-("- m+1 ) = 2" 1 " 1 = dim H Q . (37) 

The operators Bk, Ajk commute with Cj 0) ...j p _ 1 , so they leave the code subspace invariant. 
The restrictions of these operators, Bk\c and Ajk\c, satisfy the relations analogous to (j3"T|)-(|3"3D . 
Thus the correspondence Bk i— > Bk\c, Ajk ^ Ajk\c extends to a *-algebra homomorphism 
H : L(7^o) — ► L(£). But dim£ = dim7i , hence \x is an isomorphism. It can be represented as 
fJ.(X) = JX.P , where J : Ho — > C is a unitary map which is unique up to an overall phase. 

Now that the main construction is complete, we can give exact rules for converting the 
allowed (geometrically local) operations on LFMs into qubit gates. These rules are almost 
obvious but still worth putting them into a rigorous form. For most generality, consider a 
two-LFM gate application U{j, k}, where (j,k) G E. Here U is a physical operator acting 
on 2 LFMs (= 4 Majorana fermions = 2 qubits), so it can be expressed in terms of A = 
-ic c 2 = -^[Ol^tl], B' = -icQd = a 2 [0] and B" = -ic 2 c 3 = a z [l}. Applying U to LFMs j 
and k means substituting Ajk, Bj, Bk for A, B' , B" . Instead of that, we actually do another 
substitution: 

v ik : A ^ A jk , B' ^ Bj, B" ^ B k . (38) 

It extends to a *-algebra homomorphism Q — > L(B® Sjk ), where Q is the algebra generated by 
A, B' and B" (= the algebra of parity-preserving operators on two qubits), and Sj k is the set of 
qubits consisting of all edges incident to j and k. (Note that we do not have to restrict Ajk, Bj, 
B k to the subspace C because the cycle relation (^) is irrelevant in this context). So, U{j, k} 
is simulated by Vj k (U). 

It turns out that this simulation is pretty efficient even if d (= the largest degree of a vertex 
in the graph) is not a constant. W. 1. o. g. m > 2, so Sj k contains at least one edge besides 
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(j,k), say, (J, I). Then Bj, Bk and Aj k have no nontrivial relations (like BjB k = 1)- More 
exactly, the map Vj k is injective, so Bj, B k and Aj k satisfy exactly the same relations as cr z [0], 
a z [l] and — a y [0}a x [l]. It follows that there is a "symplectic transformation" of the form 
X i — > WXW^ which takes Bj, B k and Aj k to a z h a* k and —a^a^. The unitary operator W acts 
on the qubits from Sj k and can be represented as a product of elementary "symplectic gates" : 
H, A(a x ) and A(i). It is easy to show that 0(d) applications of these gates are sufficient. Hence 
executing the operator 

p jk (U) = U[(j,l),{j,k)} W (39) 
takes 0(d) one-qubit and two-qubit gate applications. 

In the above analysis, we did not take into account the number of operations required to 
prepare an initial state £ L, which is necessary to begin simulation. For definiteness, we 
will consider the unique quantum state |£) G L which satisfies additional constraints: 

B k \£) = |f) for each k. (40) 

Note that B k = —ic2kC2k+i = 1 — 2% a k, so \0 represents the fermionic state |0, . . . ,0). By a 
general argument, the qubit state |0) ® ■ ■ ■ (8> |0) can be transformed into |f ) by applying 0(u 2 ) 
symplectic gates. 

However, if we continue the speculations about fermions being possibly mimicked by bosons, 
the vacuum state of the bosonic system must absorb new degrees of freedom as the Universe 
expands. Certainly, this process should be reversible, i. e. it should also work when the space 
shrinks. In our model, shrinking the space corresponds to contracting some edges. More 
specifically, contracting an edge (j, k) means removing it while identifying the vertices j and k. 
If both j and k are connected to the same vertex I, a double edge between j = k and / appears; 
it must be then reduced to a single one. We should be able to update our "vacuum state" |f) 
through these transformations. (The qubit being removed should come out in the state |0)). 
One can show that a single contraction or reduction step requires 0(d) symplectic gates. This 
does not involve any nonlocality and can be done simultaneously in different places, which is 
quite consistent with the idea of adiabatic vacuum transformation in the expanding Universe. 



9 Quantum codes by Majorana fermions 

Some symplectic codes on qubits can be conveniently described in terms of Majorana fermions. 
We will show how it works for the Shor code |TJ]]. Whether this approach can help to find new 
codes still remains to be seen. 

By inverting the formula (H), we can represent the operators a k , a\ in terms of the Majorana 
operators: a k = \(c2k + ic2k+i), % = |( c 2fc — ic2k+i)- However, one can also introduce another 
set of annihilation and creation operators which will satisfy the same commutation relations: 

h = ^ (c T (2fc) + *C T ( 2fe+ l)) , b k = \ ( C r(2k) - ^C T ( 2 fc + l)) , (41) 

where r : {0, . . . , 2m— 1} — > {0, . . . , 2m— 1} is arbitrary permutation on 2m elements. One can 
define a quantum code by fixing the occupation numbers of some of the new LFMs, e. g. by 
requiring that each codevector \ip) satisfies b k b k \i/j) = for k = 1 . . .m — 1. In other words, the 
set of stabilizer operators is 

X k = -ic Ti 2k) c r ( 2 fc+i), fc = l,...m-l. (42) 
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The number of encoded qubits is m — (m — 1) = 1. The logical operators for this code are 
generated by c T ( ) (the encoded cr x ) and c T (i) (the encoded a y ). 

Thus each permutation r : {0, . . . , 2m— 1} — > {0, . . . , 2m— 1} defines a quantum code which 
encodes 1 qubit into m qubits. It turns out that these codes can have arbitrary large code 
distances. We are to define a family of such codes which can be considered as slightly modified 
Shor codes. Let I > 2 be an integer (the index of a code in the family), m = I 2 . There will be 
two types of stabilizer operators: 



Z k = —i C 2 H+1 C 2 {k+2)l-2 , k = . . . I — 2, 
Yk,j = —i C 2 fc/+2j+3 C2kl+2j , k = 0, . . .1 — 1, j 

They can be expressed in terms of the Pauli operators as follows : 



0,..J-2. 



(43) 



(k+2)l-2 

Z k = a x [kl]a x [(k+2)l-l] [] Y kJ = a y [kl+j]a y [kl+j+l] (44) 

s=kl+l 

For example, let us consider the I = 3 code. Its 9 qubits can be arranged into a 3 x 3 array 
as follows 

0, 1, 2, 
3, 4, 5, 
6, 7, 8. 



Then the stabilizer operators of the first type become: 



a 



I, 



a 



a 
a x 
I 



I, 



I, 

<T Z , 
O x 



(These are not matrices. We just mean that, for example, Z = o x 
I ® I ® 1). The stabilizer operators of the second type are 



Yr 



0,0 



a" 



I, 
L 



Yia = 



( 

I, 

V h 

(I, I, 

I, a y 
\I, I, 




^,1 



^2,0 — 



I, 
L 



a" 

I, 
I 



( /, 



Y 1 



( I, I, I, \ 
/, /, /, 

V <T V , <T V , I ) 



.0 



Y2A 



<J 



V I, 



I, \ 
I J 



( I, I, I, \ 
/, /, /, 



If one performs the cyclic permutation a x i— > cr y i— > a z i— > cr 11 , the operators l^j turn into certain 
stabilizer operators of the Shor code. The operators will not become exactly the same as in 
the Shor code. Each of them will contain two a y instead of a x . However, the code distance of 
this code is the same, namely, 3. 

For arbitrary /, the code distance is di = I. (The proof is essentially the same as for the 
Shor code). 
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